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Correlated-photon imaging, popularly known as ghost imaging, is a technique whereby an image 
is formed from light that has never interacted with the object. In ghost imaging experiments two 
correlated light fields are produced. One of these fields illuminates the object, and the other field is 
measured by a spatially resolving detector. In the quantum regime these correlated light fields entail 
entangled photons produced by spontaneous parametric down- conversion. To date, all correlated- 
photon ghost-imaging experiments have scanned a single-pixel detector through the field of view 
to obtain the spatial information. However, scanning leads to a poor sampling efficiency, which 
becomes worse as the number of pixels N in the image is increased. In this work we overcome such 
limitations by using a time-gated camera to record the single-photon events across the full scene. We 
produce high-contrast images in either the image plane or far-field of the down-conversion source, 
taking advantage of the EPR-like correlations in position and momentum of the photon pairs. Our 
images contain a large number of pixels, creating opportunities in low- light- level, high-resolution 
imaging and in quantum information processing. 

In a ghost imaging system the object and the imaging detector are illuminated by two correlated, spatially separated 
light beams, and the ghost image is obtained from intensity correlations between the two light fields. These optical 
correlations can occur in either the quantum or classical regime [T . Ghost imaging was first demonstrated in 1995 by 
Pittman et al. who used correlated photon pairs generated by spontaneous parametric down- conversion (SPDC). 
They used one of the down- converted fields to illuminate an aperture, and the transmission through it was detected 
by a non-imaging, single-pixel detector, termed the bucket detector. An image of the aperture was obtained from the 
coincidence counts between this bucket detector and a scanning single-pixel detector in the other beam. Bennink et 
al. [3] presented a demonstration of ghost imaging using a classical light source, and Ferri et al. [4 have since shown 
that ghost imaging using classical light sources can beat limits thought to be obtainable only when using quantum 
entangled light. Since ghost imaging was first proposed there has been a large body of published work, using both 
thermal (MTU], and quantum entangled light [6l [TIHTSl 

All previous realisations of ghost imaging using SPDC photons have relied on scanning a detector to measure the 
spatial distribution of the optical field. This scanning necessarily limits the detection efficiency of such an imaging 
system to a maximum of where N is the number of scanned pixels. In the present work we overcome this 

limitation and present a practical ghost imaging system at the single-photon level by using a time-gated, intensified 
CCD (ICCD) camera. The use of such a camera allows us to capture high- visibility ghost images across the full field 
of view. The object is placed in the one arm (object-arm) of the SPDC, the transmission through which is detected 
by a large area bucket detector. The ICCD is placed in the other arm (camera-arm) and is triggered using the output 
from the bucket detector, each frame recording the position where the single photon was detected. We compensate 
for the electronic delay associated with the triggering of the camera by means of a folded, image-preserving delay line 
as illustrated in figure [l] 

Our system was designed with interchangeable lenses that allow us to utilize either the intensity correlations of the 
photon pairs in the image plane of the down-conversion source or their anti-correlation in the far-field. This capability 
of using image-plane or far-field correlation allows us to explore the EPR-like correlations in the transverse position 
or momentum of the photon pairs [16 . 

The major advantage of full field-of-view detection, compared to a scanning system, is the dramatic increase in 
efficiency in the measurement of high-dimensional spatial entanglement. High-dimensional spatial states of single- 
photons have been used in investigations of quantum information protocols [ IT7HT9] . Recent works have demonstrated 
the potential of the multi-pixel detection of spatial entanglement [2QH25] and, albeit for multiple photons, quantum 
imaging [26 . In our work, the combination of multi-pixel detection and single-photon sensitivity across an entire field 
of view suggests applications in low-light imaging systems and quantum information protocols. 

RESULTS 

Experimental setup. A schematic diagram of our experimental setup is shown in figure [l] Our down-conversion 
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FIG. 1: Schematic diagram of experimental setup. Collinear down-converted photon pairs at 710nm are generated by 
pumping a BBO crystal with a UV laser at 355nm. Li and L3 are 50mm focal length lenses used to produce an image of 
the down- conversion source onto the beam splitter (BS). L2 is a 100mm focal length lens used to Fourier transform the down- 
conversion source onto the beam splitter. A 300mm focal length lens (not shown) is used after the beam splitter in each path 
to Fourier transform the down-converted fields at the beam splitter onto the planes of the object in one arm and onto the input 
plane of an image preserving delay line (DL) in the other arm. The image preserving delay line consists of 7 imaging systems 
with a total length of 22m. Our object (inset) is a transmissive mask with the shape of a skull printed on acetate. 



source consists of a 3mm-long beta barium borate (BBO) crystal cut for type-I phase-matching and pumped with a 
horizontally polarised, high-repetition rate laser at 355 nm. The pump laser is spatially filtered and re-collimated to 
a diameter of approximately 1.2mm (FWHM) at the crystal. The BBO crystal is oriented to provide near-collinear, 
frequency-degenerate, down-converted photons at 710nm over a lOnm bandwidth, which is determined by a high- 
transmission interference filter. We use a pellicle beam splitter (BS) to split the photon pairs into the object and 
camera arms. As we use a collinear type-I SPDC source, only 50% of the photon pairs are split into different paths 
at the beam splitter. The photons used to probe the object are detected with a bucket detector, with a field of 
view that extends over the whole scene. This bucket detector registers whether the photons have passed through the 
object whilst discarding all image information. Our bucket detector consists of lenses that couple the photons into a 
multi-mode fibre connected to a single-photon avalanche diode (SPAD). 

Our ghost image is obtained by using the bucket detector to trigger a single-photon acquisition by the ICCD 
(Andor iStar); we subsequently sum these camera images over many trigger events. To obtain in- focus images, the 
ICCD camera and the object must be located in planes where the down-converted photons exhibit strong intensity 
correlations [27^. We chose two configurations that satisfy this condition. In the first configuration, both the object 
and the ICCD are placed in image planes of the BBO crystal, where the positions of the photons are correlated. In 
the second configuration, the object and the ICCD are placed in the far-field of the BBO crystal, where the positions 
of the photons are anti-correlated, arising from their momentum anti-correlation. We choose between these two 
configurations through use of an interchangeable lens system located between the BBO crystal and the beam splitter, 
as indicated in figure [l] 

In order to allow for the electronic delay associated with the bucket detector and the trigger mechanism in the 
ICCD camera (~ 70ns), the photon in the camera-arm has to be delayed before it arrives at the ICCD. Our delay line 
(DL) consists of consecutive imaging systems with unit magnification (see supplementary figure 1). There are a total 
of 7 such imaging systems in the delay line, 4 of which are formed with 1000mm focal length lenses, and 3 with 500mm 
focal length lenses, giving a total distance of propagation of 22m. Our ghost imaging system has three key planes: 
the plane of the object, the input plane of the DL and the plane of the ICCD. These key planes are all far-field planes 
with respect to the beam splitter. Using the aforementioned interchangeable lens system, the two configurations set 
a relationship between the crystal and object/camera planes that are characterised by either an imaging system with 
magnification M = 3 or a Fourier system with effective focal length /g = 300mm (see supplementary material). 

Image acquisition. Figure |2] shows a set of ghost images taken in both the image-plane and far- field configurations 
for several different values of the delay of the triggering pulse (adjusted by changing the length of the cable between the 
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FIG. 2: Ghost images and trigger delay timing. Ghost. Images obtained in the image-plane and far- field configurations 
as a function of the trigger delay, with a total integration time of 1 hour per image. The delay is scanned in increments of 
5ns, which corresponds to adding/removing Im of cable. The images shown in (a) and (b) are enlarged versions of the images 
obtained with optimum delay time, for which the total number of photo events were 124310 and 74021, respectively. The insets 
show cross-sections of the images averaged over the 30 rows indicated by the dashed lines. 



SPAD and the ICCD). As the timing of the triggering pulse is changed relative to the time of arrival of the correlated 
photon on the ICCD the coincidence is lost and therefore no image is recorded. Figures |2]-a and[2]-b show the optimum 
ghost images in the two configurations, both with a contrast in excess of 90%. The pump power was adjusted to give 
a triggering rate from the bucket detector compatible with the maximum response rate of the ICCD. For a pump 
power of approximately 2mW the trigger rate from the bucket detector was approximately 15kHz for the image-plane 
configuration and lOkHz for the far-field configuration. All images were obtained by summing 1800 accumulations of 
2s duration. This accumulation time was chosen in order to ensure that each individual accumulation was spatially 
sparse, and therefore photon counting could be applied. To minimize the effects of noise, a photon was defined to be 
any event that exceeded both a threshold intensity and spatial extent of 2 pixels. Using this thresholding methodology, 
the average number of detected photons per frame was approximately 70 in the image-plane configuration and 40 in 
the far-field configuration. The intensity scale of the displayed images is the number of detected photons per pixel. 
No background subtraction has been applied to any of the images (see methods). 

In the image-plane configuration the intensity correlations arise from the position correlations of the down-converted 
photons, figure [2]-a. In the far- field configuration the intensity correlations arise from the momentum anti-correlations 
of the down- converted photons and hence the resulting image is inverted, figure [2]-b. As the optical magnification 
in each arm is the same, the size of both images on the ICCD match the actual size of the object. The resolution of 
both images is comparable and is set by a combination of the correlation length of the down-converted photons and 
the resolving power of our optical system (see supplementary information). 



4 



Ghost imaging 'CCD camera 




FIG. 3: Principles of a SPDC imaging system and a Klyshko advanced- wave picture. By replacing the bucket 
detector with the diode laser shown, the Klyshko advanced-wave picture enables the prediction of the spatial distribution of 
the coincidence counts that would be obtained using the SPDC imaging system. The diode laser is back-propagated through 
the system using the BBO-crystal as a partially reflecting mirror. 



Klyshko's advanced wave picture. 

By using Klyshko's advanced- wave picture [28], it is possible to predict the spatial distribution of coincidence counts 
in a SPDC experiment by using classical physics. In the advanced-wave picture, one of the detectors is replaced with a 
light source, which emits backwards through the system towards the down- conversion crystal. If the down- conversion 
crystal is viewed as a mirror, the back-propagating light is reflected at the crystal facet in the direction of the other 
detector. The overall transmission of this classical system is proportional to the coincident counts that would be 
observed in the SPDC experiment, as shown in figure [3] 

We tested our ghost imaging system according to Klyshko's picture by replacing our bucket detector with a classical 
light source, namely a laser diode at 710nm coupled to the multimode fibre. The back-propagated laser light was 
partially reflected at the crystal facet and redirected to the other arm, propagated through the delay line and recorded 
on the ICCD. We obtained a near-perfect correspondence between the back-propagated and ghost images, as can be 
seen in Figure [4] It is interesting to note that the back-propagated image in the far-field configuration is displaced 
in relation to its corresponding ghost image. This is due to the fact that the collinear phase-matching of the crystal 
results in the normal of the crystal facet differing slightly from the axis pump beam (see supplementary information) . 

An interesting feature of the back-propagated images is the fact that they are speckled. The speckle arises as a 
consequence of the first-order spatial coherence of the laser light. We have found that by shaking the fibre during 
the image acquisition time, the speckle patterns are averaged out and we obtain smooth images. In contrast, as the 
SPDC light is coherent only in second-order the ghost images exhibit a smooth intensity profile (see Figure [4|. 

Information content. The amount of information contained in the transverse distribution of a single photon 
is an important question in the context of a communication or cryptographic system [iTl [29l |30]. Two practical 
limitations to consider when using spatial states in this way are the resolution of the measurement apparatus and 
the finite size of the transverse field distribution. The resolution of our optical system in both image-plane and 
far-field configurations corresponds to approximately 20 pixels on the ICCD, equivalent to 260/im. This resolution 
was measured by acquiring ghost images of a sharp-edged object (a razor blade) and measuring the bright-to-dark 
transition of the images. The diameter of the field distributions in the image-plane and far-field plane as detected by 
the camera have a FWHM of approximately 4.3mm and 7.2mm, respectively. Taking the smaller of these distributions, 
and considering an information encoding based on binned transverse position states of the single photons, our imaging 
system can distinguish approximately 700 orthogonal states. This corresponds to a potential of 9 bits of information 
per detected photon. Figure [5] shows a set of images obtained by summing an increasing number of frames. As 
discussed above, photon-counting was performed by using intensity and spatial thresholding techniques. This gave a 
heralding efficiency of 0.2%, which was calculated as the ratio between the average number of detected photons at 
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FIG. 4: Comparison of ghost images with back-propagated images according to Klyshko's advanced wave 
picture. The speckle patterns in the back-propagated images are due to the first order spatial coherence of the illumination 
laser light, and the displacement of the far-field back-propagated images is due to the inclination of the down- conversion crystal. 
Shaking the fibre for the duration of the exposure, time averages the speckle pattern as shown. 



the camera and the triggering rate. This heralding efficiency could be increased by changing the threshold procedures 
applied, but this would come at the cost of increased dark count rates and reduced visibility We calculate the 
information content per heralded photon in our system to be log2(700 x 0.002) ^ 0.5 bits. 



DISCUSSION 



We have experimentally demonstrated a practical-ghost imaging system at the single-photon level by utilising 
a multi-pixel, time-gated, intensified camera. We recorded full-field ghost images with 90% visibility by utilising 
the spatial correlations of photon pairs in position or the anti-correlations in momentum. The ability to obtain 
high-contrast images using either position or momentum correlations is an image-based demonstration of the EPR 
phenomenon. 

We circumvent the time delay associated with the trigger mechanism of our camera and other components by 
building a free-space, image-preserving optical delay line. The ability to take full-field images greatly increases the 



6 



1 frame 
(2s) 



10 frames 




1 00 frames 


(20s) 




(200s) 



1000 frames 
(2000s) 



FIG. 5: Photon counting and heralding efficiency. Images obtained in image plane and far- field configurations for an 
increasing number of acquired frames, along with the corresponding number of detected photons. Each frame is the result of a 
thresholding procedure (see methods) of the recorded intensity distributions over an accumulation time of 2s. 



sampling efficiency. Finally, our ghost imaging system allowed us to experimentally test and validate Klyshko's 
advanced- wave picture using a full- field image. We believe this ghost imaging system could be useful both in imaging 
applications as well as in investigations of quantum information protocols with spatial states. Regarding imaging 
applications, one could envisage a non-degenerate SPDC system [31] where the object is illuminated with infrared 
light whilst the image is recorded with photons in the visible spectrum. We think this non-degenerate approach could 
be of particular interest for low-light, low-energy imaging of biological samples. 

METHODS 

Camera settings. We use an Andor iStar ICCD camera with a Gen 2 image intensifier and WR photo-cathode. The CCD head has 
1024 X 1024 pixels of size (13 x 13)/im-^. The ICCD was air cooled to —15° and operated in the direct gate trigger mode. When triggering 
in the direct gate mode, the TTL pulse from the SPAD triggers the intensifier of the camera directly, giving a gate width equal to the 
length of the trigger pulse. The maximum triggering rate of the ICCD is 50kHz when working with periodic pulses, although we work 
with a lower single count rate to minimise the possibility of overlapping trigger pulses within one response time of the intensifier. 

Threshold photon counting. Photon counting in a single-photon sensitive CCD is made by applying a threshold to each pixel 
output. By taking a set of images with the camera settings of interest but the shutter of the camera closed, we calculate the average 
dark-count rate in the camera and use this for image correction. A histogram of the pixels output from these images is used to set a 
threshold and estimate the corresponding dark-count probability. It is worth noticing that there is a trade-off between the dark-count 
rate and the detection probability; they both decrease as the threshold increases. Single-photon counting is made by thresholding each 
dark-count subtracted image. Once the threshold procedure is applied, every illuminated pixel is counted as a photon. In this regime, 
it is desirable to have a mean photon number per pixel smaller than one, since the threshold does not distinguish between one or many 
photons. 

The resolution of an ICCD is not set by the pixel size of the CCD chip, but rather by the image intensifier. This limited resolution 
manifests itself as a blooming of the photo-events across several pixels on the CCD chip. Therefore, as well as the intensity threshold, 
we spatially filter each image to remove single pixel events. Finally, in order to estimate the number of detected photons, we count the 
number of remaining events in each image. The heralding efficiency can be estimated from the average number of detected events and the 
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single count rate at the bucket detector. 
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FIG. 6: Supplementary figure. Full experimental setup of our ghost imaging system. 



Supplementary material 

Interchangeable lens systems. When exploring the position correlations of the photons we utilise lens L2 
(100mm), which performs a Fourier transform from the plane of the crystal to the plane of the beam splitter (see 
supplementary figure 1). This places both the object and the input plane of the DL in the imaging plane of the down- 
conversion source, and sets the magnification of the imaging system to \M\ = 3. When exploring the momentum 
correlations of the photons, we switch L2 for the lens system composed of lenses Li (50mm) and L3 (50mm). This 
lens system images the down-conversion source onto the beam splitter, thereby putting the reference planes in the 
far-field of the crystal. In this case, the composite Fourier system has an effective focal length of fe = 300mm. 

Image preserving delay line. We use a type-I down-conversion source with the optical axis of the BBO crystal 
orientated horizontally, thus generating vertically polarised down-converted photons. The polarisation of the photons 
in the camera-arm is rotated by 90° using a half- wave plate @ 45°, so they are transmitted through the polarising 
beam splitter (PBS). They are then imaged through the PBS to the input mirror of a telescopic imaging system with 
unit magnification constructed with Im focal length lenses. The photons propagate through this imaging system and 
are back-reflected along the same path, double passing through a quarter- wave plate set @ 45°, thus rotating the 
polarisation of the photons back to horizontal. They are therefore reflected when incident again on the PBS and are 
re-imaged on to the camera. 

Two-photon intensity correlation. In our ghost imaging system, we make use of the intensity correlations 
of the photon pairs in the image plane and far- field of the down-conversion source. The correlation length of the 
photon pairs in these planes can be estimated from the two-photon wave function and from our optical system. In 
the thin-crystal, monochromatic and paraxial approximations, the post-selected two-photon wave function of SPDC 
can be written as [32l |33] 



where qi (q2) is the transverse component of the wave vector of photon 1 (2). The function V is the angular spectrum 
of the pump beam, which is transferred to the two-photon wave function [32 , and 7 is the phase-matching function 
of the SPDC. Since the angular spectrum of the pump beam V is much narrower than the phase matching function 
7, the photon pairs exhibit anti-correlation between their transverse momenta. The two-photon wave function in 
transverse position representation is given by the Fourier transform of Eq. (IT]) : 



where pi and p2 are the transverse coordinates in the detection planes of photons 1 and 2, W is the transverse field 
distribution of the pump beam and F is the Fourier transform of the phase matching function. As the function F 
is much narrower than W, the transverse positions of the photons are correlated. We define the correlation lengths 



^(qi, q2) = v(qi + q2)7(qi - q2), 



(1) 



^(pi, P2) = w(pi + P2)r(pi - P2), 



(2) 
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of the photon pairs as the FWHM of the two-photon detection probabihties in image-plane and far-field, which are 
proportional to the widths of |rp and |Vp respectively. Using the parameters of our SPDC and optical system, we 
estimate these correlation lengths to be 



which are both below the resolution of our optical system of 260/im. 

Far-field back-propagation image offset The offset of the far-field back-propagated image on the ICCD is 
approximately 440 pixels or 5.72mm. Using the effective focal length /g = 300mm of the Fourier system, we calculate 
the angle between the normal of the crystal facet and the pump beam to he = 5.72mm/ ^ 1.1°, which corresponds 
to the inclination angle measured with the rotation stage. The inclination of the crystal does not give a displaced 
image in the image-plane configuration as the crystal is imaged to the ICCD. 




62/im, 
78/im, 
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